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I. INTRODUCTION 


The Advanced Helicopter Electromagnetics (AHE) Industrial Associates Pro- 
gram continues its research on variety of main topics identified and recom- 
mended by the Advisory Task Force of the program. The research activities 
center on issues that advance technology related to helicopter electromagnet- 
ics. While most of the topics are a continuation of previous works, special 
effort has been focused on some of the areas due to recommendations from 
the last annual conference held at Arizona State University on October 23- 
24, 1991 and the biannual meeting of 1992 held at the Naval Air Warfare 
Center on May 13-14, 1992. The main topics addressed in this report are: 

• Composite Materials 

• Antenna Technology 

The area of Composite Materials continues getting special attention in 
this period. The research has focused on: 

a. Measurements of the electrical properties of low-conductivity materials 

b. Modeling of material discontinuity and their effects on the scattering 
patterns (RCS) 

c. Preliminary analysis on interaction of electromagnetic fields with multi- 
layered graphite fiberglass plates 

d. FDTD modeling of fields penetration through composite panels of a 
helicopter. 

In previous research, our work on RCS and radiation pattern predictions 
was primarily limited to high conductivity materials which make little dif- 
ference from perfect electric conductors. Considering that low conductivity 
materials are frequently used in modern aircraft industry, we have inves- 
tigated the effects of those low conductivity materials on RCS of complex 
structures including conducting and composite materials. The spectral do- 
main Green’s function approach has been employed to model the composite 
by surface impedance, which in turn is coupled into the NEC code for pat- 
tern predictions. Some of the results are reported in this period. This report 
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also provides new measurement results on GDS composite material that we 
believe to be more accurate. 

In light of the anisotropic property of a graphite fiberglass composite, 
basic analysis has been performed on the interaction of EM waves with an 
anisotropic multi-layered composite. The analysis is aimed to lay the basis 
for the introduction of anisotropic surface impedance. From the prelim- 
inary study, it is found that low frequency fields may penetrate through 
multi-layered composites depending on the number of layers and the fiber 
orientations. 

As a follow-up, the shielding problems in a helicopter were examined 
further by using FD-TD method when some of the helicopter panels were 
replaced with composite materials. 

In the area of Antenna Technology, significant progress has been made 
in investigations of different antenna elements such as stacked microstrips, 
cavity-backed microstrips, varactor-tuned microstrip patches and ferrite- 
loaded cavity-backed slots. These novel microstrip antenna concepts show 
some attractive features overcoming certain intrinsic short-comings of con- 
ventional microstrip antennas. 

In particular, a two-dimensional analytical model of a magnetically- 
tuneable cavity-backed slot antenna is introduced in this report with specific 
intent of analyzing the scattering characteristics and the tuning capabilities 
of such an antenna at VHF and UHF frequencies. This analysis is hoped to 
demonstrate the tangibility of such an antenna and to bring a substantial 
amount of insight into the specific design considerations. 
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II 


COMPOSITE MATERIALS 


A. Introduction 

The progress in the area of Composite Materials is significant and diverse. 
The major achievements are represented by the RCS prediction of low con- 
ductive composites, also EM wave penetration through composite panels into 
a helicopter, and analysis of EM wave interaction with anisotropic composite 
plates. 

Since low conductivity materials are used frequently in modern helicopter 
platforms, it is necessary to study the effects of these materials on RCS of 
complex structures that include several types of conducting, semi-conducting 
and nonconducting materials. In this report, we investigate the RCS patterns 
of low conductivity composite materials including modeling the discontinu- 
ities between these materials and PEC. RCS comparisons were initially made 
because the effects of materials and material discontinuities are usually more 
evident in scattering patterns than they are in antenna patterns. 

To illustrate the capabilities and potentials of the FDTD method to elec- 
tromagnetic penetration and interaction problems on complex airframes, the 
perfectly conducting NASA helicopter airframe is replaced by a helicopter 
structure that is made of conducting and composite materials. The compos- 
ite materials (fiberglass and teflon) used to replace some of the helicopter 
panels are typical samples of materials being used in the construction of 
modern helicopter airframes. The helicopter is then exposed to a plane wave 
and the electromagnetic fields penetrating the structure were computed us- 
ing the FDTD method. The solution of the problem reveals the insight into 
wave penetration and interaction with on-board structures where vulnerable 
electronic equipment is housed. 

Furthermore, in this report we also present an analytical formulation 
and numerical program for systematic prediction of EM wave reflection from 
and transmission through composite plates involving stratified anisotropic 
graphite-epoxy fiber sheets. Various numerical simulation results of the re- 
flection and the transmission characteristics of stratified composites versus 
frequency are reported. A significant penetration of an electromagnetic wave 
through a composite plate at a low frequency range has been revealed by the 
numerical results. The strong dependency of the reflection and the trans- 
mission characteristics of a composite plate on the orientation of the fibers 
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composing the plate can be observed. 


B« RCS of Low Conductivity Composite Materials 

1. Introduction 

In our previous reports we presented theory and measurements of: 

a. Antenna and RCS patterns on composite materials of high conductivity 
such as graphite epoxy and fiberglass with screens near the surface. 

b. Antenna patterns on composite materials of low conductivity including 
microwave absorbers. 

In the above cases, the surface impedance of the composite material is 
evaluated using a spectral domain Green’s functions approach and utilized 
to model the structure under investigation as a wire grid. The pattern of 
the structure is then computed using the wire grid option of NEC code. 
Monopole antennas on both high and low conductivity materials have been 
built and tested. Good agreement between predictions and measurements of 
radiation patterns of these antennas have been observed. 

So far our work on RCS of composite materials has been limited to high 
conductivity materials. The RCS patterns of structures made of high conduc- 
tivity materials were very similar to those made of perfect electric conductors 
(PEC). In this report, we investigate the RCS patterns of low conductivity 
composite materials including modeling the discontinuities between compos- 
ites and PECs. 


2. Theory and measurements 

It has been shown in our previous reports that the surface impedance Z 5 of 
a two-layer medium including a layer of composite material (Figure la), can 
be written as 


(i) 

Gs = T lj Z u 


Z a = ( Te + ZtGs)(Yt + TjG,)~ l 
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where Yt, Te , Tj are the transmission matrices of the composite 
material and the subscript L denotes the medium below the composite ma- 
terial. At the limit when the thickness d of the bottom layer equal zero, the 
surface impedance becomes 


^ ^ ^ — i 

Z s = ZtTj (2) 

To verify our theory, a 4.0’ , x2.7”x0.06” aluminum plate was partially cov- 
ered by a layer microwave absorber as shown in Figure 2 and the RCS pattern 
of this structure was measured at 4 GHz (the lowest operating frequency of 
our anechoic chamber at the time of the measurements). The dimensions 
of the absorber were 2.0”x2.7”x0.03” (half the size of the plate) The elec- 
tric properties of the absorber material were measured at X band (8.2-12.4 
GHz) using the waveguide technique and are shown in Figures 3-6. The 
measured properties were then interpolated to e r =14.9-j0.1, // r =2.0-jl.6 at 
4 GHz (RCS measurement frequency). There has not been any other ac- 
curate technique, available to us, to directly measure the absorber material 
properties at 4 GHz. The surface impedance was computed using the above 
theory for different angles of incidence. Very small variations of the surface 
impedance versus the angle of incidence have been predicted (41.8 - j48.2fi 
±10%). Using this computed value of the surface impedance, the NEC code 
was used to model and compute the RCS pattern of the structure shown in 
Figure 2. The results are shown in Figure 7. Note that, since the aluminum 
plate in this was relatively thick (0.06”), it was supposed to model the plate 
by two grid layers at the surface of the plate. However, to save CPU time 
and computer memory, only one PEC grid at the center was used to model 
the aluminum plate. The measured pattern and the calculated pattern using 
our theory are in good agreement. The absorber has a noticeable effect on 
pattern symmetry as well as on the intensity of the main beams. The main 
beam on the absorber side was about 3 dB less than the PEC side. 

3. Conclusion and future work 

The effects of low conductivity composite materials on the RCS patterns 
of antennas have been analyzed and presented. The code has been vali- 
dated by computing the surface impedance comparing predicted patterns of 
a plate partially covered with a microwave absorber material to measure- 
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ments. Presently, we are investigating the scattering from more complex 
structures such as helicopters and other complicated platforms. 

C. Electromagnetic Penetration through Composite 
Panels 

1. Simulating EMI problems of the helicopter scale model with 
the FD-TD method 

To simulate electromagnetic interaction and penetration problems, some of 
the perfectly conducting panels of the helicopter are replaced by composite 
material panels. Specifically, at the front section of the helicopter, teflon with 
e r = 2.10 - J0.014 is used. At the bottom of the helicopter, four panels are 
replaced by fiberglass with t r = 4.60 - J0.017 at 500 MHz. The configuration 
of the helicopter used in this simulation is shown in Figure 8, where the 
composite material is shown as the shaded region. The composite materials 
used in this case are representative of materials used to construct modern 
helicopter airframes. At the frequency of 500 MHz, the composite material 
panels are assumed to have very small thickness compared to the wavelength 
of operation (typically A/500) and so they are modeled as thin plates. 

A plane wave traveling parallel to the length of the helicopter (positive 
x-axis), with vertical polarization ( E z and H y components) is incident on 
the hybrid (perfectly conducting and composite) helicopter airframe. The 
amplitude of the incident electric field is assumed unity. The plane wave 
incident on the helicopter structure is scattered by the helicopter in all dif- 
ferent directions. Because of the presence of some composite material panels, 
a significant amount of energy penetrates through the composite panels. To 
study the amount of penetration through the airframe, a horizontal cut was 
taken at a z = constant plane of the helicopter model shown in Figure 8. The 
cross section of the horizontal cut taken is shown in Figure 9. The helicopter 
wings are shown as grid meshes. The incident plane wave striking the heli- 
copter airframe is also shown in the figure. Over the horizontal cut shown 
in Figure 9, the amplitude distribution of the penetrating fields can be ob- 
tained. The FDTD solver allows the option to specify horizontal and vertical 
cuts in the problem space and obtain the amplitude and phase distribution 
of preselected field components. In this case, the magnetic field amplitude 
distributions for H x and H y are obtained from the FDTD solver. 


6 


In the previous progress report the same problem was solved, however, 
with the direction of incidence of the plane wave being in the negative x- 
axis. In this progress report the incidence angle is changed so that the plane 
wave hits the structure from the opposite direction. As will be demonstrated 
shortly, the penetration fields are smaller than in the previous case. 

The horizontal cut taken is further divided into three separate regions (A), 
(B), (C) and the fields in each region are plotted separately. The results of the 
y- and x-components of the magnetic field distribution inside the helicopter 
are illustrated in Figures 10 and 11, respectively. These are normalized with 
respect to the amplitude of the incident magnetic field and are plotted in 
dB (20 * log l0 {H total / H inc )) over the helicopter cross section. Because the 
incident magnetic field component is polarized in the y direction, the H y 
field (co-polarization) distribution is much higher than the corresponding 
H x field (cross-polarization) distribution. The fields penetrating through the 
helicopter in this case are smaller than the case considered in the previous 
progress report, where the plane wave struck head-on the composite material 
section of the helicopter structure. 

2. Conclusions and future work 

Some of the perfectly conducting panels on the helicopter structure are re- 
placed by composite materials and the interaction of the composite airframe 
with a plane wave has been examined. The replacement of the perfectly con- 
ducting panels with composite ones allows the penetration of the plane wave 
inside the helicopter structure which acts like a cavity. The distribution of 
penetration fields has been examined by obtaining the x- and y-components 
of the magnetic-field distribution inside the helicopter. 

Future work will include more elaborate study of electromagnetic pene- 
tration into more complex helicopter structures by considering inside wires, 
metal supports, electronic systems, or electronic housing departments. Ways 
to obtain the current distribution from the magnetic field distribution on the 
helicopter airframe are also under consideration. 

In addition, radiation patterns on the helicopter structure for frequen- 
cies lower than 200 MHz, for example 20 MHz, will be attempted using the 
FDTD. In low frequency applications the grid size selection in the FDTD 
modeling has to be extremely small (less than A/100). This small grid size 
creates problems in the FDTD modeling. One problem is that in FDTD sim- 
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ulations a distance of about one A is left between the structure being modeled 
and the outer boundary of the FDTD problem space where absorbing bound- 
ary conditions are applied. This requirement will have an enormous increase 
in the memory needed to solve the radiation problems of low frequencies. 

One way to reduce the one A distance requirement is by using more ac- 
curate absorbing boundary conditions. The accuracy of these conditions 
depends on the order of approximation used. At present the second-order ap- 
proximations were considered to be sufficiently accurate for the applications 
that were attempted. In the next reporting period progress on increasing 
the order of absorbing boundary conditions of up to seventh order for three 
dimensional problems will be reported. Hopefully, these high-order bound- 
ary conditions will allow the distance between the structure and the outer 
boundary to be reduced to 0.1 A, for example, and thus allow modeling of low 
frequency radiation problems. 

D. Interaction of EM Waves With a Stratified Com- 
posite Plate 

1. Introduction 

The interaction of electromagnetic waves with fiber-reinforced matrix com- 
posites has gained more attention of the research community in both ma- 
terial and electrical science [1, 2]. The considerable interest has also been 
manifested in research areas of astronautics, electromagnetic engineering and 
space vehicle design [3, 4, 5]. In fact, there are many situations where the 
consideration of this interaction becomes inevitable, such as the radiation 
effect of composites on helicopter antennas, and electromagnetic shielding of 
equipment on an airplane by composite fuselage. 

More specifically, graphite-epoxy composite plates are widely used to 
build helicopters and fighter airplanes due to their superior mechanical 
strength and light weight. Composite materials are anisotropic and less con- 
ductive than a metal. Therefore, when a helicopter is illuminated by an EM 
wave, it is easier for the EM wave to penetrate into the interior of the heli- 
copter to cause EM interference problem. This interference may disturb the 
on-board electronic and computer systems. A more serious problem of the 
same nature is encountered when the helicopter is hit by lightning. 

The complexity of a basic analysis is simply attributed to the anisotropic 
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electromagnetic property of the composites. Physically, the anisotropy is di- 
rectly related to the orientation of fibers composing the composites. Thus the 
three principal axes of the permittivity and permeability, which characterize 
the electromagnetic properties of the composites, are found to be in paral- 
lel or perpendicular to the fiber orientation. The analysis in this report is 
based on the fundamental assumption that the permittivity and permeability 
tensors of a composite are representable by diagonal matrices. 

The goal of this preliminary research is to study the electromagnetic prop- 
erties of composite materials over a wide frequency band, and understand the 
basic electromagnetic properties of composite materials by analyzing the EM 
scattering behavior, or radar cross section, and the EM shielding capability 
of anisotropic composites. In the next research period, a proper surface im- 
pedance model is to be established, which will be coupled into the NEC and 
ESP codes by using the impedance loading concept. 

2. Theoretical model 

To simplify the analysis, we first consider a composite laminate constructed of 
M plies, each ply consisting of unidirectional fibers embedded in an organic 
matrix. The fiber orientation with respect to the x-axis is denoted by an 
angle 0, as shown in Figure 12. The composite is illuminated by a linearly 
polarized plane wave at normal incidence. The polarization angle is S as 
shown in Figure 12. 

In an anisotropic ( conductivity cr ) medium, the Maxwell’s equations are 
represented as 


V x E = -jutiH (3) 

V x H = cr E + juj e • E ; = j u t ■ E (4) 

where the dyadic permittivity is defined as 

- a 
£ = - — + £ 

e is the dielectric permittivity, and /< is the permeability. Taking the curl of 
(1) yields 

V xV x E = “jw/i xH 
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By using (4) and vector identity, we have 


V (V • E) - V 2 E = -ju(j u 7 -E) 

It can be shown that in the composite 

V • E = 0 

Thus, the wave equation in the anisotropic medium takes the form of 

V 2 E + u> 2 p 7 ■ E = 0 (5) 

In the principal coordinate system of the dyadic permittivity, the permit- 
tivity can be written as a diagonal matrix 


e = 


e; 0 0 

0 e* 0 
0 0 c! 


where z-direction is defined as the propagation direction of the incident wave, 
p and q are the other two principal directions, and 


e; = + e 

p ju 


s = 




+ c 


c* = -4 + € 

In the ( p, q ) coordinate system, (3) can be decomposed into two components: 

V 2 E P + to 2 p e p Ep = 0 (6) 

V 2 E q + c o 2 pt q E q = 0 (7) 

Since the incident wave is a plane wave propagating in a normal direction to 
the composite, the above equations simplify to 


<P 2 

— Ep + u) ptp Ep = 0 


(8) 


10 



( 9 ) 


(P 9 

~ Eg + V 2 (it q Eg = 0 

az z 

The homogeneous solution to (8) and (9) is provided by 

E p = E+ e~ jkpZ + E~ e ]hpZ (10) 

E q = E+ e~ jkqZ + E~ e jkqZ (11) 

where the superscript + denotes the forward propagating wave and the — 
denotes the reflected wave, and k p and k q are wavenumbers in p-axis and 
< 7 -axis directions, which are related to electric parameters of the medium in 
which the wave is propagating by 

k 2 p = u 2 fit p (12) 


k 2 = u 2 (ie q (13) 

Once the solution of the electric fields is obtained, the magnetic fields are 
deduced from 

II = — V x E 

Ufi 


Thus, we obtain 


If we define 


H p = ^-(E+e~ jkqt - E~e jkqZ ) 

Ufl q q 

(14) 

H q = -2L (£+ e ~ ikpZ - E~ e jfcpZ ) 

UJ jj 

(15) 

H v = H* e~ jkqZ + H; e jkqZ 

(16) 

H q = H+ e~ jkpz + H~ e jkpZ 

(17) 


then, the amplitudes of the magnetic fields of the reflected and transmitted 
waves can be related to those of the electric fields by a matrix of 


K1 

1 


0 

1 ^ 

H t J 

_ 

0 

\f^P _ 

U + 


( 18 ) 
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and 


’ b; ' 

1 

l^-J 

v // 7 



' E-' 


[e; J 


(19) 


0 

0 

Consider a plane wave incident on the interface between two media as 
shown in Figure 12. The components of the electric fields in the medium a 
with respect to the global coordinate system ( Xi, X 2 ) can be projected on 
to the principal coordinate system of the medium a by 


’ E\' 


l^J 



co$0 a — sin9 a 
sin6 a cos9 a 

and the magnetic fields are projected as 


' B\ ' 


. H 2 . 



cost) a 
sin0 a 


— sin9 a 
cos9 a 


E\ 


Hi 


H] 


<? j 


by using (18), (21) reduces to 


B\' 


COS 9 a 

— sin9 a 

1 

Hi. 


sin9 a 

COS 6 a 

s/B 


& 

0 


Hi 


cos9 a 

— sin9 a 


sin9 a 

cos 9 a 



0 



'Ey 


.E}\ 


-E\- 


E p J 


' H\ ' 

1 


^ . 


cosO a — sind a 
siri0 a cosB a 


0 

L\/5 0 


cos9 a sin9 a 


' El' 

— sin9 a cos9 a 


E\ 


Hi 

m 


Na 


E\ 


(20) 


( 21 ) 


( 22 ) 


(23) 


(24) 


(25) 


where N a is called the characteristic matrix for the medium a. Following the 
same procedure, similar relations for the reflected and transmitted waves can 
be established, that is, 


m 


= - N a 


El 

El 


(26) 
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(27) 


H[ 

m 


N b 


E{ 

El] 


with Nt is the characteristic matrix of the medium b. 

Without the loss of generality, the interface of the media a and b is as- 
sumed to be located at the plane of z = 0 . Matching the tangential com- 
ponents of the electric and magnetic fields on the interface leads to 


and 


E\ 

E\ 


H\ 

m 



- E[ - 


' El' 

■ 

[El 


[El J 

1 

' HI ‘ 


’ Hi' 


. m _ 


H l J 


Using the above results, we can write it out 

E\ 


N a 


E j 
Ei 


- N b 


El 


= N a 


Solving simultaneously (28) and (30) leads to 


El 

El 


and 


El 

El 


N a + N b ]- l [N a - N b ] 


E \ 
El 


E\ 

El 


= 2 [N a + N b ]~ l N a 


E\ 


(28) 

(29) 

(30) 

(31) 

(32) 


In the preceding discussions, the wave is incident from one side only. As 
shown in Figure 13, if waves are incident from both sides of the interface, 
their reflected and transmitted counterparts are given by 


E+ *= E+ — E r + (33) 

E- = Ef - E r - (34) 

It is helpful to define the effective reflected and transmitted waves here as 


P + = E+ + E; (35) 

P- = E t “ + E r + (36) 
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Since the reflected waves can be characterized by reflection coefficients as 


E+ = r • Ef (37) 

E; = - r e; (38) 

where r is the reflection dyadic from the medium a to the medium 6, which 
can be represented by the characteristic matrix as 

r= [N a + N b ]- l [N a - N b ) 

Consequently, the effective reflected and transmitted electric fields can be 
represented in terms of the reflection dyadic by 


P+ = Ef + ( — r * E~ ) = Ef + E+- r E” 


) • Ef - r ■ E" 

(39) 

Er+ ? E,+) 

-E- - E r "+ r -E t + 

(40) 

( / - 7 ) • E- + ? • E+ 

(41) 


Up to this point, we have considered only the case of reflection by and 
transmission through one interface between two media. If the composite 
is a multi-layered structure, the solution of the electric fields in different 
media can be obtained in a similar way. Considering two adjacent interfaces 
as shown in Figure 14, we use the subscript m — 1 to denote the interface 
between the (m — 1 ) th ply and mth ply. Inside the mth ply of the composite, 
the electric fields at the (m — 1 ) th interface are related to that at the mth 
interface by a transmission matrix A m as 


where 


A 


m 


(E, + )m = A m P+ 
(Ej ) m ~ A m P m 

cos9 a — sin9 a 
sinO a cos 0 a 


^jkpdrn 

0 


0 

—jkqdtn 


(42) 

(43) 


cos0 a sin0 a 
— sin9 a cos9 a 


(44) 
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If there are M layers in the composite, based on the previous formulations, 
the effective reflected and transmitted electric fields inside the mth ply can 
be represented by that inside the two adjacent plies as 

P+ = ( 7 + r m ) A m —\ P+_, - A m P- +1 (45) 

and 

P m = ( I ~~ r m ) A m P^_|_1 + r m ^4m-lP m _i (46) 

Since the excitation wave is incident on to the first interface, applying 45 
and 46 at the first interface yields 

Pf = (7 + ?i) E inc - A x P 2 - (47) 

Pr = (7 - r x )A x p-+ ?1 E inc (48) 

or 

Pf - Pr + ^lP 2 " = Einc (49) 

Pf + r x A x PJ = (7 + ?i) • Einc (50) 

where E mc denotes the excitation electric field. Similarly, at the interface of 
m = M + 1, we obtain 

P M+1 = (7 + r M+l ) ^mPm (51) 

PtfH — r M+\ AmPXt (52) 

or 

AmP\ { — Pjf + i + Pjif + i = 0 (53) 

r M + 1 -4a/ Pm — p m+i = 0 (54) 

In general, at the mth interface, we have 

(7 + r m )/l m _ 1 P+_ ] - P+- r m ^ m p- +1 = 0 (55) 

r m — P m + (I + r m )v4 m P m+1 = 0 (56) 

In fact, the effective reflected and transmitted electric fields also can be 
represented by a cascaded matrix form that relates the fields at one interface 
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to that at another interface. Thus, the solution of the total reflected and 
transmitted fields can be achieved by inverting the cascaded matrix. 

Consider the ( M + l)th interface, from equations (51) and (52), we have 


P M - [(/ + r M+l)^Af] ’Pm+I 


and 


p m+ 1 = r M+i A m [(I + r m+i) Am] lp M+i 
which can be written in a matrix form as 


= Zm + i 


p+ 

r M + 1 


(57) 

(58) 


, 

T) 


[( / + r \f+ 1 Am] 1 

0 


p+ 
r M+ 1 

[ Pm+i 


r m + i Am[{I + r 

0 


0 


(59) 


Equation (59) is a special form at the (M -f l)th interface. In general, 
the effective reflected and transmitted electric fields at one interface are rep- 
resented by that at the next interface as 


P+ . 

m-l 

P" 


= z, 


M 


p £ + 1 

P Af— 1 


(60) 


where 


with 


a 

b 

c 

d 


Z\f = 


a b 
c d 


[( / + r mMm — l] 1 

[( / T r m ) A m _i] A m 

^ m — I T 771 ) A 772 — 1 j 

(/ ^ m ) -Tm ^ m — 1 [( I T m ) m — 1 ] ^ m 


Based on these equations, we can relate the total reflected and transmit- 
ted electric fields to the external excitation as 


E 

J - # me 


— i 

+ 

i 

+ 

►— * 

i 

p t 

p r . 


. r, (7 + r,)- 1 (7 - ro-M,+ (7 + ?,)-> 7, ^ _ 

. P 2 . 


( 61 ) 
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(62) 


E- 

J - J me 

pr 


— Zq ■ Z\ ■ Zi 


Zm+i 


P+ 

r M+l 


’ E mc ' 

— 7 . 

P + 

. p r . 

— £J 

0 


E 



*21 


p+ 

r M+ 1 

. p r . 


. Z2X 

*22 _ 


0 


Thus, the total reflected and transmitted electric fields are given by 


(63) 
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and 

Pjr+i = z n * 'E'inc (66) 

The total reflectance R and transmittance T are then determined by 


R = 
T = 


fP'r I 2 

|E mc p 

(67) 

|Pt/ + il 2 
i e,„ c i 2 

(68) 


Equations (59) to (68) are fundamental relations used to analyze the 
reflection and transmission characteristic of a stratified composite illuminated 
by plane electromagnetic wave at normal incidence. Extensive numerical 
simulations have already been conducted. Typical examples are illustrated 
in the following section. 


3. Numerical examples 

This section presents various numerical simulations of the interaction of a 
plane wave with a multi-layered composite structure. The electric parameters 
are cited from reference [1], and are assumed to be valid in the frequency band 
of interest. 

First, we try to obtain the scattering characteristics of a one-layered 
prepreg AS/3501-6. Different polarizations of the incident electric fields 
with respect to the fiber running direction are used to calculate the total re- 
flectance and transmittance. Figure 15 shows the reflectance of a one-layered 
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prepreg before curing when an incident electromagnetic wave is illuminated 
at relative polarization angles of 0°, 45°, and 90°. As the frequency of the 
incident wave is increased, we observe that the reflectance at 0° approaches 
unity very rapidly. However, the reflectance at the other polarization angles 
is relatively small, even when the frequency is up to 4 GHz. Thus, we see 
that the reflectance is strongly dependent on the polarization of the incident 
field. Figure 16 shows the transmittance of the specified structure. At 0° 
excitation, the transmittance is very small when the frequency is above 1 
GHz. At 90° excitation, the transmittance is much larger. Especially at low 
frequency end, it is very close to unity, which implies the transparency of the 
structure to the incident wave. 

Figures 17 and 18 show the results of two-layered structure that consists 
of two plies of prepregs placed so that the fiber orientations of each ply are 
in parallel. The same polarization-dependency of the incident field is also 
observed. The penetration of the electric field, which is perpendicular to the 
fiber orientation, is also dominant especially at the low frequency range. 

Figures 19 and 20 show the comparison of the re flectances and transmit- 
tances between different layers of composite prepregs. The fiber direction of 
one prepreg is perpendicular to that of the adjacent prepregs. The curves are 
shown for two-layered, eight-layered and sixteen-layered structures which are 
illuminated by a 45°-polarized electromagnetic wave. As we observed, when 
the frequency is beyond 0.6 GHz, the transmittance is very close to zero. In 
other words, if the frequency is higher than 0.6 GHz, the wave penetration 
is negligible. But, at the low frequency range, the penetration can be over- 
whelming even though sixteen-layered plates are used. For a helicopter to 
prevent the penetration of lightning, special construction is needed. 

Figures 21 and 22 show the reflectance and transmittance of two-layered 
and 45° crossed prepregs when the incident electric field is polarized at dif- 
ferent angles. It is seen that once there is a ply that has a fiber direction in 
parallel to fhe polarization of the incident electric field, the transmittance of 
the entire plate becomes very small at frequencies higher than 1 GHz. The 
strong dependency on the polarization of incident fields is again observed. 
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Figure 1: Modeling geometry of a composite plate 
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Figure 2: Modeling geometry of a PEC partially covered by GDS 
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Figure 3: Measured real part of the permittivity of GDS 
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Figure 4: Measured imaginary part of the permittivity of GDS 
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Figure 5: Measured real part of the permeability of GDS 


23 


9.6 10.3 11. 11.7 12.4 

FREQUENCY (GHz) 





FR 







Figure 8: Partially-conducting/partially-composite material structure used 
to study electromagnetic penetration into the helicopter airframe. 
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Figure 9: Helicopter cross section at a constant z plane. The magnetic field 
distributions in regions (A), (B), (C) are obtained from the FDTD solver. 
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Figure 10: H v magnetic field distribution inside the helicopter structure at 
500 MHz. 
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Figure 11: H x magnetic field distribution inside the helicopter structure at 

500 MHz. 
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Figure 12: Incident, reflected and transmitted fields by an interface between 
media a and b. 
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Figure 13: Effective reflected and transmitted fields on an interface between 
two media. 
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Figure 14: Reflection and transmission of the electric fields by multi com- 
posite layers 
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Transmittance 



Figure 16: The transmittance of a one-ply AS/3501-6 composite plate illu- 
minated by an incident wave polarized at 0°, 45° and 90° 
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Figure 18: The transmittance of a two-ply AS/3501-6 composite plate illu- 
minated by an incident wave polarized at 0°, 45° and 90° 




Reflectance 



Figure 19: The reflectance of multi-layered perpendicularly crossed AS/ 350 1 - 
6 composite plate illuminated by an incident wave polarized at 45° with 
respect to the fiber orientation of the prepregs 
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Figure 20: The transmittance of multi-layered perpendicularly crossed 

AS/3501-6 composite plate illuminated by an incident wave polarized at 45° 
with respect to. the fiber orientation of the prepregs 
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Reflectance 



Figure 21: The reflectance of a two-plv 45° crossed AS/3501-6 composite 
plate illuminated by an incident wave polarized at 45° with respect to. the 
fiber orientation of the face prepreg 


Transmittance 



Figure 22: The transmittance of a two-ply 45° crossed AS/3501-6 composite 
plate illuminated by an incident wave polarized at 45° with respect to. the 
fiber orientation of the face prepreg 
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Ill- ANTENNA TECHNOLOGY 

A. Introduction 


During the past quarter, significant progress in four areas of conformal an- 
tenna technology has been made. These areas are: 

a. stacked microstrip patch antennas, 

b. cavity-backed microstrip patch antennas, 

c. varactor-tuned microstrip patch antennas, and 

d. ferrite-loaded cavity-backed slot antennas. 

Both microstrip and cavity-backed slot antennas are widely used at UHF and 
above in aerospace applications such as helicopters because they are confor- 
mal to the vehicle surface as well as light weight, low cost and easy to man- 
ufacture. The first three topics mentioned above represent novel microstrip 
antenna concepts that have the potential to overcome certain intrinsic disad- 
vantages of microstrip antennas such as narrow bandwidth and the tradeoff 
between scan volume and bandwidth in large arrays without significantly 
compromising the attractive features that make microstrip antennas so de- 
sirable. In this section, progress to date in each of these areas is discussed, 
and plans for future work are outlined. 

B. Stacked Microstrip Patch Antennas 

It is well-known that a tradeoff exists between bandwidth and scan volume 
in a microstrip patch element phased array. The most straightforward way to 
improve the bandwidth of a microstrip patch antenna (for a given substrate 
dielectric constant) is to increase the separation between the patch and the 
ground plane. However, an unloaded grounded dielectric substrate supports 
tightly bound surface wave modes which become more strongly excited as 
the substrate becomes thicker. In a microstrip antenna, these surface wave 
modes do not contribute to main beam radiation and hence represent a loss 
mechanism. In a large phased array of microstrip patch elements, guided 
w T ave effects related to the surface wave modes of the unloaded substrate can 
lead to large mismatches (scan blindnesses) at certain scan angles. Thus, 
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while use of a thicker substrate results in an improvement in bandwddth, the 
variation in active impedance versus scan becomes more rapid. Two ways to 
overcome this tradeoff between scan volume and bandwidth in large arrays 
of microstrip elements have recently been proposed. The first is to stack 
the patches in two dielectric layers [6j. By stacking patch antennas, one 
can improve the impedance bandwidth of the antenna and/or provide for an 
antenna which can operate at two widely separated frequencies. The second 
is to break the substrate into discrete regions by placing metal walls between 
the patch elements [7]. The generic geometry of a stacked microstrip patch 
antenna is shown in Figure 23. Computer programs for the following cases 
of infinite arrays of probe-fed stacked patch antennas have been successfully 
developed and validated: 

a. one-probe-fed rectangular patches 

b. two-probe-fed rectangular patches 

c. one-probe-fed circular patches 

d. two-probe-fed circular patches 

As an example of the results generated by our computer codes for in- 
finite arrays of stacked patches, consider Figure 24. This figure illustrates 
broadside-matched reflection coefficient versus scan angle in both E- and 
H-planes for a one-probe-fed stacked rectangular patch infinite array. The 
polarization mismatch for the array is shown in Figure 25. The computer 
code required 35.2 hours of CPU time to run on a VAX 3100 workstation. 

During the next quarter, additional calculated results will be generated 
with particular attention paid to the two-probe-fed geometries for polarization- 
diverse microstrip antennas. We will also begin to extend this work to the 
case of single stacked patch elements. 

C. Cavity-Backed Microstrip Patch Antennas 

Another way to overcome the tradeoff that exists between bandwidth and 
scan volume in a microstrip patch element phased array is to surround each 
patch element with a cavity. Doing so breaks the substrate into discrete 
regions, and prevents the excitation of surface wave modes. Recently on 
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a project for the Engineering Foundation, we investigated a cavity-backed 
patch geometry in which probe-fed circular patches are backed by rectangular 
cavities [7], [8]. It was demonstrated that a dramatic improvement in the scan 
performance of microstrip antenna phased array printed on a given substrate 
can be achieved by the use of a metallized cavity-backing. 

It has recently been suggested that the use of circular cavities may be 
more appropriate due to the ease with which they may be machined [9]. In 
addition, the use of alternate feeding techniques such as aperture coupling 
is highly desirable. During the past quarter, an investigation of probe-fed 
circular patches backed by circular cavities was begun. We have successfully 
derived the dyadic Green’s functi on fo r the circular cavity in the form of 
a discrete vector Hankel transform (DVHT) which should be the basis for 
development of an extremely efficient numerical code. During the next quar- 
ter, we will begin to develop numerical codes for this geometry by modifying 
existing codes for probe-fed circular patches backed by rectangular cavities. 

D. Varactor- Tuned Microstrip Patch Antennas 

Narrow bandwidth is a major disadvantage of conventional microstrip and 
cavity-backed slot antennas. This disadvantage precludes their use in a num- 
ber of applications that could benefit from antennas possessing their many 
desirable characteristics. In this and the next section, we discuss some of the 
efforts underway to develop tunable conformal antennas. The basic concept 
is to start with an inherently narrow-bandwidth antenna and somehow vary 
one of its parameters allowing the antenna to be tuned over a fairly wide 
bandwidth. In many applications, it is actually more desirable to use such 
a tunable, inherently narrow-bandwidth antenna than an inherently wide- 
bandwidth antenna because of the ability of the narrow-bandwidth antenna 
to reject out-of-band interference. 

In past quarterly reports, the basic operation of the varactor- tuned mi- 
crostrip patch was explained, theoretical results for the input impedance 
and radiation pattern of circular patches fabricated with both silicon and 
GaAs varactor diodes were presented, and experimental results for the input 
impedance of circular patches fabricated with silicon varactor diodes were 
presented. During the past quarter, experimental results using GaAS var- 
actor diodes were obtained. Figure 26 shows both measured and calculated 
input return loss versus frequency as a function of bias level for the GaAs 
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varactor-tuned patch. 

Using varactor diode tuning, we have successfully designed frequency- 
agile microstrip patch antennas. While both silicon and GaAs diodes enable 
tuning of the patch frequency response, our calculations reveal that GaAs 
diodes are a better choice because the higher Q of the diode dramatically 
improves the efficiency of the antenna. Calculated minimum return loss fre- 
quencies and the antenna efficiencies (due to diode loss) at these frequencies 
as a function of bias level are summarized in the following table. 


Varactor 

Bias Level 

(V) 

min. R. L. freq. 
(GHz) 

efficiency 

(%) 


no diode 

2.82 

100 

Silicon 

20 

2.66 

89.6 


4 

2.18 

22.5 


no diode 

2.82 

100 

GaAs 

20 

2.74 

99.3 


4 

2.52 

83.2 


2 

2.28 

55.3 


A paper [10] summarizing our work to date with varactor-tuned mi- 
crostrip patches is to be presented at the upcoming IEEE Antenna and 
Propagation Society International Symposium in Chicago. 
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E. Ferrite- Controlled Cavity-Backed Slot Radiators 

1. Introduction 

The concept of using ferrite to magnetically tune cavity-backed slot (CBS) 
antennas has been around for quite some time. Previous literature, however, 
fails to bring to light any substantial experimental work which verifies the 
capabilities of such an antenna. Moreover, to the best of the authors knowl- 
edge, no analytical work has been presented which attempts to model the 
characteristics of a ferrite-tuned cavity radiator. 

During the past quarter, a two-dimensional analytical model of a magneti- 
cally-tuneable cavity-backed slot antenna was developed with the specific 
intent of analyzing the scattering characteristics and the tuning capabilities 
of such an antenna at VHF and UIIF frequencies. The work presented here, 
discusses the tangibility of such an antenna and brings a substantial amount 
of insight into the specific design qualities needed to realize these goals. 

2. Background 

Ferrite substrates are ferrimagnetic materials constructed from solid ceramic 
materials which have been sintered with certain metal oxides at high tem- 
peratures. At microwave frequencies they exhibit strong magnetic effects, 
which result in anisotropic behavior. By applying an external DC magnetic 
field, the permeability tensor of the ferrite is altered, and thus the electrical 
characteristics of the ferrite substrates used in the CBS structure are varied. 
Because of this, the antenna will resonate at different frequencies for different 
values of the applied DC magnetic field. 

There are two types of modes that can be excited inside the cavity. The 
first type of mode is the dynamic mode, which depends mainly on the di- 
electric properties of the materials inside the cavity, including dielectric con- 
stants of those materials and their relative permeabilities. The second type 
of mode is the magnetostatic mode. For magnetostatic modes, the perme- 
ability tensor gives rise to low frequency magnetostatic volume-wave modes, 
which become appreciable with large external biases. It is these volume wave 
modes that allow a microwave CBS antenna to operate in a frequency range 
lower than those of the dynamic mode. 
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3. Two-dimensional analytical model 

The analysis presented in this report utilizes a spectral domain admittance 
Green’s function in conjunction with the transmission matrix method to ac- 
count for multiple substrate layers within the cavity. In this method, closed 
form expressions of the Green’s function are derived for the regions above 
and within the cavity, and take into account the anisotropic characteristics of 
the ferrite. The full-wave formulation uses the physical equivalence principle, 
along with Galerkin’s procedure, to obtain the scattered fields produced by 
the antenna when the structure is excited by an incident plane wave. Figure 
27 shows the schematic diagram of a multiple-layered cavity filled with both 
dielectric and ferrite materials. Tuning capabilities of both magnetostatic 
and dynamic modes were analyzed for a variety of layering configurations 
and different cavity dimensions. The results presented here examine several 
antenna characteristics which may have an appreciable effect on the perfor- 
mance of the antenna and its tuning capabilities. 

4. Dielectric-loaded cavities vs. ferrite-filled cavities 

In order to design a physically realizable magnetically-tunable CBS antenna, 
the primary question that needs to be addressed is: What effect does the 
ferrite have on the cavity-backing which will allow for tuning capabilities, 
and what characteristics will allow it to operate at lower frequencies? In 
order to answer this question, we examine two cases. The first case is a 
2D rectangular cavity-backed slot in an infinite ground plane partially filled 
with a dielectric substrate. The second case is developed by replacing the 
dielectric with a ferrite substrate of the same relative dielectric constant. 

If we illuminate the CBS structure by an incident plane wave with hard 
polarization (TMz), we obtain the resulting scattering width of the antenna 
as a function of frequency given in Figure 28. The peaks in the figure repre- 
sent the modes at which the antenna will resonate. 

It is evident from the figure that by replacing the dielectric with a fer- 
rite substrate, the resonant frequency of the dominant dynamic mode of the 
structure has been shifted as a result of the alteration of the electrical di- 
mensions of the cavity. Even more interesting, is the presence of a weak 
magnetostatic volume-wave mode at approximately 60 MHz for H o —0 Oer- 
sted. The presence of this mode coincides with the findings of the classic 


46 



papers by others on the propagation of volume- wave magnetostatic modes in 
ferrite substrates. We can see from the figure, hovrever, that without an ex- 
ternal bias, this mode is far too weak to be considered as a useable resonant 
mode. 

In order to utilize these volume-wave modes as a source of radiation, it is 
necessary to apply a substantial DC magnetic field to the ferrite substrates. 
By increasing the magnetic bias, the magnetostatic modes become gradually 
stronger in magnitude and shift over in frequency as well. Figure 28 shows the 
results of the scattering width of the antenna as the external bias is increased 
from 0.0 Oersted to 300.0 Oersted. As shown in the graph, the volume-wave 
modes tune over a substantial frequency range whereas the dynamic mode is 
less affected by the magnetic bias. The following table shows the magnitudes 
and resonant frequencies of the magnetostatic volume-wave modes as the 
external bias is increased. 


DC Bias (Oersted) 

Frequency (GHZ) 

Scattering Width (dB/m) 

0.0000 

0.2007 

-26.73 

100.00 

0.7570 

-8.56 

300.00 

1.399 

-5.42 


5. Design considerations 

5. a. Airgap and low dielectric material effects In previous work on 
ferrite antennas, there were no air gaps between radiating elements (such as 
slots, microstrip patches, etc.) and ferrite layers. Also, the ferrite has always 
been metal backed, i.e. there was no separation between the ferrite substrate 
and the ground plane. Clearly, losses in ferrite will have a significant impact 
on the total antenna efficiency if there is no separation between the radiating 
elements and the ferrite. Also, ground plane has the effect of attenuating 
volume wave modes and therefore, there has to be a separation between 
ferrite and ground plane to support low loss volume waves. 


47 



First, we consider the effects of increasing the airgap between the slot and 
ferrite in order to increase the antenna efficiency. To investigate such effects, 
the thickness of the ferrite and the low dielectric substrate at the bottom 
of the cavity were assumed constant and the airgap is varied. The analysis 
shows that the scattering width of the volume-wave modes increases as the 
air gap increases, as shown in Figure 29. Increasing the air gap also resulted 
in increasing the bandwidth of both the magnetostatic and dynamic modes 
of the cavity. 

Secondly, we study the effect of placing a low dielectric substrate at the 
bottom of the cavity to shield the ferrite from the ground plane. From 
Figure 30, it is evident that an increase in the ground layer will increase the 
scattering width of the magnetostatic modes of the antenna by approximately 
1.0 dB. From these results, we conclude that the inclusion of a sufficient 
airgap and a low dielectric grounding layer will add approximately 4 dB to 
the scattering width of volume- wave modes. 

5.b. Effect of cavity width It is well known that the cavity dimensions 
will have a drastic effect on the dynamic resonant modes of the CBS antenna. 
It is not well known, however, what effect these dimensions will have on the 
magnetostatic resonant modes of the ferrite-loaded cavity antenna. From the 
2D analytical model, we can show that an increase in the cavity dimensions 
along the direction of magnetization will in fact greatly increase the magni- 
tude of the ferrite modes as well as shift them tow’ard lower frequencies of 
resonance. 

Figure 31 shows the results of the scattering width of the antenna as the 
cavity width is varied from 0.6 inches to 2.0 inches. As the figure shows, 
an increase in the cavity width has a large impact on the magnitude of 
the magnetostatic mode but not on the dynamic mode of the cavity. The 
following table shows the actual increases in the magnitudes of the lowest 
order volume- wave mode and the associated frequency shift. 
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Cavity Width (inches) 

Frequency (GHZ) 

Scattering Width (dB/m) 

0.6 

0.8013 

-11.88 

1.0 

0.7520 

-8.25 

2.0 

0.6286 

-4.78 


5.c. Effect of the slot width The width of the slot is yet another cavity 
dimension that can be considered in t he design of a CBS antenna. An increase 
in the slot dimension may greatly increase the bandwidth of the dynamic 
resonant mode, but efficiency is reduced as a result. The magnitude of the 
magnetostatic mode, on the other hand, is not significantly altered, but does 
cause a slight shift in the frequency of resonance, as shown in Figure 32. This 
shift will only become appreciable for large cavity dimensions, where the slot 
width is allowed to vary over a significant range. 

F. Future Work 

We are in the process of building an optimized magnetostatic volume wave 
antenna. We acquired strong permanent magnet tips to provide the required 
magnetization for the ferrite layers. Other materials, needed to build the 
antenna, are already available to us. We intend to extend our analysis to 
include the three dimensional case and evaluate the difference between the 
two dimensional and the three dimensional theories. 
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Figure 23; Geometry of a stacked patch phased array. The bottom patch is 
probe-fed via the center conductor of a coaxial connector, and the top patch 
is parasitic. 
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Figure 24: Broadside-matched reflection coefficient magnitude versus scan 
angle for an infinite array of stacked one-probe-fed rectangular microstrip 
patches. The substrate parameters of this array are: 

t r = 2.5, d\ = 0.0510, t T = 1.5, d 2 ~ 0.0510, where A 0 is the free-space 
wavelength. The lower patch is 0.29A o by 0.29 Aq ; the upper patch is 0.285A o 
by 0.285A o . The lower patch is fed by a probe feed of radius 0.01450Ao at a 
point 0.0725A o from the center along the x-axis. The element spacing in the 
x-direction is 0.51A o ; the element spacing in the y-direction is 0.5A o . 
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Figure 26: Measured and calculated input return loss of the GaAs varactor 
tuned microstrip patch antenna for bias voltages of 20, 4 and 2 volts. 





Figure 27: Schematic diagram of a multiple-layered cavity filled with both 
dielectric and ferrite materials. 
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Figure 28: Scattering width of a CBS antenna versus frequency. 
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Figure 29: Scattering width of the volume- wave modes versus frequency and 
air gap. 
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Figure 30 : Scattering width of a CBS antenna with ground layer. 
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Figure 31: Scattering width of a CBS antenna versus the cavity width and 
frequency. 
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Figure 32 : Scattering width of a CBS antenna versus frequency and the slot 
width. 
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